Crystallization and phase-separation in non-additive binary hard-sphere mixtures 
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We calculate for the first time the full phase-diagram of an asymmetric non-additive hard-sphere 
mixture. The non-additivity strongly affects the crystallization and the fluid-fluid phase-separation. 
The global topology of the phase-diagram is controlled by an effective size-ratio j/cff, while the 
fluid-solid coexistence scales with the depth of the effective potential well. 
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Entropically driven phase transitions have received 
much attention lately because of their direct relevance to 
the observed phase-behaviour of colloidal suspensions [Q . 
Theoretical work has focussed on the asymmetric binary 
hard sphere (HS) system, a deceptively simple mixture 
of large and small particles, which exhibits an interesting 
competition between demixing into dilute and concen- 
trated suspensions of large particles, driven by the famil- 
iar osmotic depletion effect Q, and freezing into an or- 
dered crystalline phase. Recent Monte Carlo (MC) simu- 
lations 1^ of binary colloidal dispersions of spherical par- 
ticles with additive diameters ui and (72, and size-ratio 
y = oijox < 0.2, have convincingly demonstrated that 
the demixing transition conjectured earlier [Q is always 
preempted by a direct freezing of a low concentration dis- 
ordered ( "fluid" ) phase into an FCC crystal of large par- 
ticles. It was argued elsewhere § that a small degree of 
non-additivity of the diameter ai2 , determining the dis- 
tance of closest approach between large and small parti- 
cles, might drive the demixing transition from metastable 
to stable. This view has been supported by recent MC 
simulations |^ which show that non-additivity signifi- 
cantly lowers the packing fraction at the critical point 
of the fluid-fluid coexistence curve. However, to reach 
firm conclusions concerning the observability of a fluid- 
fluid demixing transition, the effect of non-additivity on 
the freezing transition must also be considered explicitly, 
in order to map out a complete phase diagram. This 
is precisely the objective of this letter, where the im- 
portant influence of non-additivity on the global phase 
behaviour of highly asymmetric hard sphere mixtures 
is evaluated within a systematic Statistical Mechanical 
perturbation treatment. Even a small non-additivity is 
shown to have a large effect on the interpretation of the 
phase behaviour, and experiments on sterically or charge- 
stabilized binary "HS" colloids are shown to generically 
exhibit non-additive behaviour. 

Consider a binary system of HS with distances of clos- 
est approach (Tq/3 (1 < a,/3 < 2), such that cth = ci, 
^22 = ^2, and: 

CT12 = ^((Tii +(T22)(1 + A), (1) 

where the non-additivity parameter A can be positive or 



negative. The case A = corresponds to additive HS, 
and has been widely studied with the usual techniques 
of the Statistical Mechanics of fluids. An extreme exam- 
ple of non-additivity is provided by the case (T22 ~ 0, 
A > 0, which is the familiar Asakura-Oosawa model 
of colloid-polymer mixtures, allowing for full inter- 
penetrability of polymer coils. This system is predicted 
to exhibit phase separation into dilute and concentrated 
colloid fluid phases for large enough A |0,D, in good 
agreement with experimental findings The present 
work is concerned with the case of small positive non- 
additivity for HS mixtures with small size-ratios which 
would lead to a metastable demixing transition in the ad- 
ditive limit, A = 0. A > will obviously favour demixing 
since phases involving a majority of particles of the same 
species will allow a more efficient packing . 

To examine the possibility of a fluid-fluid phase separa- 
tion for A > 0, it seems natural to develop a perturbation 
theory expansion in powers of A around a reference mix- 
ture of identical composition, involving additive HS. The 
free energy per unit volume /o = Fq/V, of the latter is 
rather accurately given by the semi-empirical equation of 
state proposed by Mansoori et. al. pT| ], which improves 
on the compressibility equation of state derived from the 
analytic solution of the PY equations Note that 

neither predicts a spinodal instability for any size-ratio 
or composition, as characterized by a vanishing determi- 
nant of the stability matrix M = \d^ fQ/duadUf^l, where 
ria = Na/V is the number density of a spheres. How- 
ever for fixed packing fraction, they predict that M ap- 
proaches zero as suggesting that for larger and 
larger size-asymmetry, smaller and smaller perturbations 
to the purely additive case can drive the stability matrix 
negative. The additive case is in a sense "marginal" to 
phase-separation, which explains why the very existence 
of such a separation is so sensitive to the particular ap- 
proximations used 0. 

A direct application of standard thermodynamic per- 
turbation theory to the unlike pair potential: 

^i2''('') — '''i2'(''/(l + '^A)), where A = and 1 corre- 
spond, respectively, to the additive reference system with 
diameter af^^ , and to the non-additive mixture of inter- 
est, leads to the first order correction to the free energy: 
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In Eq. ^, {af2 ) is the contact value of the unlike pair 
distribution function of the reference mixture taken from 
the analytic solution of the PY equation for an additive 
binary HS mixture p2[. 
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FIG. 1. The upper curves are fluid-fluid spinodals for 
y = 0.2, A = 0.033 as obtained from the 4 theories de- 
scribed in the text: Mansoori reference (o), PY reference 
io), Barboy-Gelbart {+), and BPGG (*). The lower 3 
sets of curves compare the fluid-solid binodals calculated 
by 1st order perturbation theory in the potential of Eq. 
(^) (dashed lines) with MC results (solid lines) for the 
additive case, A = 0. From top to bottom, y — 0.2, 0.1, 
and 0.05 respectively. 

The fluid-fluid spinodal curves resulting from this two- 
component perturbation theory using both the Mansoori 
and PY jl^ reference Fq's are compared in Fig. |l| 
to an approximate series expansion due to Barboy and 
Gelbart |Jl|l for a size-ratio y = 0.2 and A = 0.033. 
Note that the results are plotted for a semi-grand ensem- 
ble: the binary mixture is in osmotic equilibrium with a 
reservoir of small spheres which fixes their chemical po- 
tential fi2'i the thermodynamic variables controlling the 
phase behaviour are the packing fractions, rji = 7rnicrj'/6, 
of the large spheres in the mixture and 773 of small spheres 
in the reservoir (or equivalently (12)- In sharp contrast 
to the additive case where different theories yield dissim- 
ilar results 1^, all 3 approaches yield similar results for 
the spinodal curves. Moreover, the are consistent with 
results from numerical solutions of the BPGG integral 
equation at least for low values of rji. As expected, 
the small correction due to the non-additivity triggers a 
phase-separation which is absent in the two additive ref- 
erence systems (A = 0). The demixing transition, which 
is marginal for additive HS's [Q, is strongly enhanced 
by a modest degree of non-additivity. Moreover, increas- 
ing A shifts the demixing transition to lower and lower 
packing fraction r/2, of small spheres, as shown in Fig. 

The above two-component perturbation scheme cannot 
be adapted to investigate the freezing of non-additive HS 



mixtures, mainly because the crystal phase of the ad- 
ditive reference system is poorly understood for small 
size-ratios y. HS alloys form substitutionally disordered 
crystals for y > 0.85 Il6|, and interesting super-lattice 
structures for y < 0.6 [^|, but the structure for smaller 
values of y (say y < 0.3) is not well understood, although 
it has been conjectured that the large spheres might form 
an FCC lattice, permeated by a fluid of small spheres, 
at least for sufficiently small y Q. To avoid these dif- 
ficulties, one may resort to an effective one-component 
description, by integrating out the degrees of freedom 
of the small spheres for any given configuration of large 
spheres. This procedure leads to effective interaction po- 
tentials between the large spheres which are determined 
by the free energy (or grand potential) of the inhomo- 
geneous fluid of small spheres in the "external field" of 
the larger particles. The basic philosophy behind effec- 
tive potentials in complex fluids is that the (considerable) 
initial effort spent deriving accurate effective potentials 
is recouped when they are used as input into the well de- 
veloped machinery of liquid state theory . In the case 
of colloid-polymer mixtures, this approach leads to the 
Asakura-Oosawa effective pair potential between large 
spheres |19 , and for additive HS mixtures the effective 
pair potential has recently been extended to include ex- 
cluded volume correlations between small spheres [^^D . 
The procedure used in these references may be gener- 
alised to derive the following effective pair potential be- 
tween large spheres in a non-additive HS mixture: 

PVcs{r) =00 ; r < 0-1 

+ {'n-2? [10 h{rf- 12 yh{r)]\ ; ai < r < ai(l + yeff), 



where the effective size-ratio is: 
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while the function h{r) — (1 + j/off) ~ f/<^i- The term 
linear in 7/2 is the purely attractive Derjaguin form of the 
Asakura-Oosawa potential, with the effective size-ratio 
y^g, while the higher order terms describe the partially 
repulsive effects of the correlation-induced layering of the 
small spheres around the large spheres. As shown in the 
inset of Fig. ||, increasing A at fixed 773 deepens the at- 
tractive well, while the correlation-induced repulsive bar- 
rier remains roughly the same. 

The effective interaction for r > cri is treated as a per- 
turbation of the one-component HS reference system, and 
standard first order perturbation theory JT^,]?! is applied 
to calculate the free energy and determine the fiuid-solid 
phase coexistence. The predictions of this effective one- 
component perturbation theory are compared in Fig. ^ 
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to MC data for the additive (A = 0) case (which in turn 
compares well with full 2-component MC simulations [^); 
the agreement is seen to be good, thus lending confidence 
to the extension of the theory to the non-additive case, 
for which no MC data is available. Interestingly this suc- 
cess is not mirrored for the fluid-fluid transition: the 1- 
component perturbation theory predicts a critical point 
at much higher values of 771 than was found in the sim- 
ulations, demonstrating that care must be taken when 
applying ideas culled from the field of simple fluids to 
the effective potentials arising from complex fluids. On 
the other hand, the two-component perturbation theory 
embodied in Eq. turns out to be much more accurate 
and has been used throughout to determine the demixing 
of non-additive HS's into two fluid phases. 

Resuhs for 0.02 < A < 0.25 are shown in Fig. |. 
Both the fluid-fluid spinodal and the liquidus line of the 
freezing transition are seen to shift to lower 772 as A in- 
creases, so that the fluid-fluid coexistence curve remains 
metastable up to A i:2 0.2, at which stage the critical 
point becomes stable. The two-component perturbation 
theory cannot be trusted for significantly larger values of 
A, because of its first order nature. On the other hand 
non-additivity is expected to remain rather small for ster- 
ically or electrostatically stabilized colloids, so that the 
present calculations imply that it may be unlikely to ob- 
serve fluid-fluid phase coexistence in asymmetric colloidal 
mixtures. 
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FIG. 2. Fluid-fluid spinodals (dashed lines) and 
fluid-solid binodals (solid lines) for y = 0.2, and 
A = 0.02,0.033, 0.05, 0.1, &0. 25 (from top to bottom). 
The long-dashed curve is the fluid-solid binodal for 
A = 0. Inset: effective potentials of Eq. (^) for 
y = 0.2,7)5 = 0.3, A = (solid line) and A = 0.05 
(dashed line) 

However, the calculations lead to another unexpected 
prediction. If the phase diagrams in Fig. ^ are rescaled, 
by choosing the reduced well depth at contact /3e = 
/?ycff(cri), rather than 772, as the thermodynamic vari- 
able along the y axis, the fluid-solid liquidus curves for 
all values of A fall practically on top of each other. This 
"quasi-universal" behaviour is shown in Fig. 0: the fluid- 



solid liquidus line scales with well-depth at contact, and 
is largely independent of the shape or range (character- 
ized by ?/off) of the effective pair potential between the 
larger spheres. 
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FIG. 3. Same as Fig. ^ but with the y-axis re-scaled 
w.r.t. the well-depth at contact of Eq. (H). Inset: two 
typical effective potentials from Eq. dsl) for yoff ~ 0.2, 
and ?/cff = 0.44; both result in very similar liquidus lines! 

This behaviour contrasts sharply with the fluid-fluid 
coexistence curves, which are very sensitive to the range 
and shape of the effective potential. Note that the fluid- 
fluid demixing transition becomes stable for j/off ^ 0.4, 
which is comparable to the results of direct simulations 
of the Asakura-Oosawa model for colloid-polymer mix- 
tures |l|- 

The possibility of phase-separation in binary hard- 
sphere mixtures inspired a series of experiments on hard- 
sphere like-colloids The hard-sphere repulsion is 
achieved by steric stabilization, with a repulsive outer 
brush of co-polymer, or else the colloids are stabilized 
by repulsive electrostatic interactions, with enough salt 
added to obtain a very short screening length. Interac- 
tions between charge-stabilized binary colloids are well 
described by the DLVO potential j2|j2j]. For sterically 
stabilized colloids, the interaction between two spheres 
can be estimated by combining the interaction between 
2 flat brush-covered surfaces, which is well described 
by the Alexander-de Gennes theory |2^], with the Dcr- 
jaguin approximation, which results in: V2i{t — 022) = 
7/yii(r-crii), and yi2(r-cri2) = {1y j {\^y))Vi_\{r- axx)-, 
where Vap{r — (Jap)) is the interaction potential between 
the brushes of species a and (3. The repulsive interaction 
is smaller between smaller particles because the brush 
surface overlap area is smaller as a function of distance 
r. The 3 effective hard-sphere diameters can be estimated 
by the diameter at which the interaction is a few times 
fcsT, and then used to determine the non-additivity. For 
both sterically and electrostatically stabilized binary col- 
loidal mixtures this results in: 
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where / is the decay length of the brush |23| , or the Debye 
screening length of the DLVO interaction. To first order, 
the non-additivity depends neither on the pre-factors in 
the interactions, nor on exactly how many times fcsT 
is chosen as the hard-sphere criterion; only the "soft- 
ness" of the potential enters. Since the logarithmic fac- 
tor in Eq. (||) is negative, sterically or electrostatically 
stabilized colloidal mixtures will manifest negative non- 
additivity. For example, the popular model colloids made 
of PMMA (polymethylmethylacrylate) cores are stabi- 
lized by a PHS (poly-12-hydroxysteric acid) brush, typi- 
cally 10 — 12 nm wide, with / w 3 — 4 nm ||2^, a number 
confirmed by direct measurements with a surface-force 
apparatus [pSj . Thus a mixture with the larger species at 
400 nm and the smaller at 40 nm (y^O.l) will have a neg- 
ative non-additivity of about A = —0.01; i.e., the cross- 
diameter (712 is 1% less than the additive cross-diameter 
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FIG. 4. Fluid-solid binodals for y = 0.1075, A = 
(solid lines) and A = —0.01 (dash-dotted lines) plotted 
for absolute 772, (not r;^) to compare to the liquidus line 
from Imhof and Dhont for the same y (dashed line). 
Fig. ^ demonstrates that even a small amount of neg- 
ative non-additivity has a relatively large effect on the 
phase-diagram, raising the liquidus line to higher effec- 
tive values of 772 . No fluid-fluid de-mixing is expected for 
A < 0. 

In conclusion then, we have demonstrated that non- 
additivity has an important effect on the phase-diagram 
of asymmetric binary hard-sphere mixtures, affecting 
both the fluid-fluid and the fluid-solid lines. The fluid- 
solid liquidus lines show a near-universal behaviour, and 
experimental hard-sphere like colloidal dispersions are 
expected never to separate into two fluid phases. 
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